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In non- cent rosy mmetric magnets, the chiral Dzyaloshinskii-Moriya exchange stabilizes Skyrmion- 
strings as excitations which may condense into multiply modulated phases. Such extended 
Skyrmionic textures are determined by the stability of the localized solitonic Skyrmion cores and 
their geometrical incompatibility which frustrates regular space- filling. We present numerically exact 
solutions for Skyrmion lattices and formulate basic properties of the Skyrmionic states. 
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INTRODUCTION 

During the last years the investigation of chiral mod- 
ulations and Skyrmionic states, both in bulk non- 
centrosymmetric magnetic materials and nanomagnetic 
systems (see e.g. [3, 0, 0, 0, 0, E| and bibliography 
in 0, 0]), has become a hot topic. Recently different 
specific modulated states have been discovered in non- 
centrosymmetric magnets and multiferroics including the 
manifestation of Skyrmionic states in the cubic helimag- 
net MnSi [1, 0, Q. Such Skyrmionic states have been 
predicted two decades ago [8|, and subsequent theoret- 
ical studies within the standard Dzyaloshinskii theory 
Q have revealed a number of novel modulated states 
in magnetic materials with intrinsic and induced chi- 
rality [3, 0, E3, El- In this contribution, we outline 
main results on chiral Skyrmionic states in the low tem- 
perature region [To[ and near the ordering temperature 
|2| derived within the standard phenomenological the- 
ory Q . We also present new numerically exact solutions 
for Skyrmion lattices. We formulate basic properties of 
Skyrmionic states and elucidate physical mechanisms un- 
derlying their formation and stability. 

CHIRAL FLUX-LINES: THE BUILDING BLOCKS 
OF SKYRMIONIC MATTER 

According to Dzyaloshinskii [9], the magnetic energy 
density of a non-centrosymmetric ferromagnet can be 
written as 

w = A ^(diMj) 2 - M • H + w D (M) + w (M), (1) 

and includes the exchange stiffness with constant A, the 
Zeeman energy, and chiral Dzyaloshinskii-Moriya (DM) 
coupling (wd)- (M) combines internal energy contri- 
butions independent on gradients of the magnetization. 

Chiral modulations arise as a result of the competition 
between the exchange and DM interactions [9[ , and a suf- 
ficiently strong magnetic field compresses them into local- 
ized two-dimensional vortex-like structures (Fig. [T] (a)) 



[1, • These chiral axisymmetric structures have a non- 
trivial topology and smooth non-singular cores with defi- 
nite sizes. The DM couplings are described in continuum 
models of magnetic materials by so-called Lifshitz invari- 
ants, energy contributions linear in first spatial deriva- 
tives of the magnetization, A-^ = Mid^Mj — Mjd^Mi 
[9|. For MnSi and other B20 compounds belonging to 
the crystallographic class T, wd = D (A^ + A^J + A^) 
=D M-rotM (where D is a Dzyaloshinskii constant). For 
uniaxial non-centrosymmetric classes the functionals wd 
are given in Ref. 

For the sake of simplicity we include into model 
(PQ) only these basic interactions essential to stabilize 
Skyrmionic states, and neglect other less important en- 
ergy contributions (as magnetic anisotropy, stray-fields, 
magneto-elastic coupling). 

For low temperatures M = const, and the equations 
minimizing energy (pQ) with wo = include solutions for 
axisymmetric localized structures ijj = ip ((/)), = 0(p), 
where M = M(sin cos -0; sin sin ip\ cos 0)) and the spa- 
tial variable r = (pcos </?; psin(/?; z) is written in cylindri- 
cal coordinates (Fig. [T]). The solutions _0(</>) are deter- 
mined by crystal classes of the system [8[. Particularly, 
ip = (j) for C nv symmetry, ip = <p — tt/2 for D n and cubic 
T classes, and ip = tt/2 — <j) for D2d symmetry (Fig. HJ 
(a-d)). The polar angle 0(p) is derived from equation 

, fd 2 IdO 1 . . \ 
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sin 2 sin = 0, (2) 

p 2M ' u 

with the boundary conditions 0(0) = 7r, 0(oo) = 0. This 
yields the solutions for isolated Skyrmions which are ra- 
dially stable at well-defined sizes L oc \D\/H [8[ (Fig. [1] 
(f,g)). The chiral interactions play the crucial role to sta- 
bilize Skyrmions. In centrosymmetric systems (D = 0) 
such solutions are radially unstable and collapse sponta- 
neously under the influence of the applied magnetic field 
or intrinsic short-range interactions 

Commonly in nonlinear field models Skyrmions arise 
due to specific invariants described by higher order 
spatial derivatives (so-called Skyrme mechanism). In 
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FIG. 1: Cross-section through an isolated Skyrmion in ferromagnets with Cnv symmetry (a). Projections of the magnetization 
onto the basal plane in the Skyrmion core for ferromagnets with C nv (b), D n (c), and D^d (d) symmetries. A hexagonal 
Skyrmion lattice for ferromagnets with C nv symmetry (e). Evolution of the magnetization 6(p/R c ) with increasing magnetic 
field (f) (adopted from Ref. [13]). Equilibrium size R c of the Skyrmion cell as a function of an applied magnetic field 
(Ld = 4:7V A/ D) (g). Inset shows the distribution of the magnetization in the Skyrmion lattice for H/Hd = (I), 0.6 (II), 0.8 
(HI) (H D = MD 2 /(2A)). 



condensed-matter systems there are no physical interac- 
tions providing such energy contributions. Hence, chi- 
ral couplings present the unique mechanism to stabi- 
lize Skyrmionic textures in ordered condensed matter 
systems described by a large class of nonlinear field 
models |2[. This singles out chiral condensed-matter 
systems with Lifshitz-type of invariants (including non- 
centrosymmetric magnets, multiferroics, and chiral liq- 
uid crystals) into a particular class of materials with 
Skyrmionic states. The Skyrmions in chiral magnets can 
be thought of as isolated filaments within spatially ho- 
mogeneous phases. Isolated Skyrmions remind Abrikosov 
vortices in type II superconductors or thread-like textures 
in nematic liquid crystals [2[ . Contrary to these defected 
patterns with singularity in the core, the distribution of 
the order parameter in Skyrmions is smooth (Fig. [T]). 



SKYRMION LATTICES VERSUS HELICOIDS 

For strong DM interactions isolated Skyrmions con- 
dense into lattices (Figs. □ (e), [2] (c-f)). These 2D 
modulated textures are alternatives to common one- 
dimensional (helical) modulations. In the model with a 
fixed magnetization modulus appropriate for low temper- 
ature the equilibrium parameters of the Skyrmion lattice 
have been calculated in circular cell approximation 
At zero field the equilibrium period of the lattice is close 
to the helicoid period (2R C « Ld)- Near the critical 
field H c = 0.8132 the lattice transforms into a system of 



isolated Skyrmions by infinite growth of the period and 
localization of the core (Fig. ffl(f), (g)) [Tol ] . 

Near the ordering temperature, the magnetization 
modulus becomes small (M <C 1) and strongly de- 
pends on the applied field and temperature. A brute- 
force minimization of the functional (pQ) with wq(M) = 
a(T — T C )M 2 + bM A in zero field yields solutions for a 
"half- Skyrmion" staggered structure and hexagonal lat- 
tice with radial variation of the magnetization modu- 
lus in the lattice cells ( Fig. [2]). [2]. The Skyrmion 
lattices are characterized by a strong variation of the 
cell sizes and transformation of their structures near cell 
boundaries. However, they preserve axisymmetric distri- 
bution of the magnetization near the cell center. This 
remarkable property is due to specific energetics of the 
Skyrmions. " Double- twist" rotation of the magnetization 
near the Skyrmion core leads to larger energy reduction 
than in " single- twisted" helical phases while edge areas of 
the cell have larger energy density than the helical states 
[H . This explains the unusual axial symmetry of the cell 
cores and their stability. The condensation of Skyrmions 
creates spatially inhomogeneous twisted phases as a re- 
sult of space-filling by these multi-dimensional solitonic 
objects. 



SUMMARY 

Skyrmionic textures reveal common features imposed 
by the fundamental physical mechanisms underlying 
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FIG. 2: (Color online) A "half-Skyrmion" lattice cell (a). Profiles 6i(p) and M(p) (b) in a hexagonal cell for two different 
directions i through the core (dashed-dotted line (i — 1) and dotted line (i — 2)) are plotted together with the corresponding 
profiles 9o(p) and Mo(p) for the circular cell approximation (solid line). Inset in (b) shows Vi = \0% — 0o|/#o as functions of p. 
Contour plots of M z (x,y) (c,d) and M(x,y) (e,f) for a square (c,e) and hexagonal (d,f) lattices. The results are derived for 
reduced parameters of model JT]) near the ordering temperature: a = a{T - T C )A/D 2 = 0.23, b = bA/D 2 = 0.05. 



their stability. Summarizing the findings of previous pa- 
pers d, H, [Io[ and results of this contribution we for- 
mulate basic properties of extended states built from 
chiral Skyrmion solutions, (i) In magnetic materials 
lacking inversion symmetry chiral DM interactions sta- 
bilize isolated Skyrmions, axisymmetric localized struc- 
tures with a fixed rotation sense and definite shape and 
sizes, (ii) Strong DM coupling leads to the condensa- 
tion of Skyrmions into lattices (Figs, dj E])- The effec- 
tive interactions between the Skyrmions are weak rela- 
tive to energies determining their radial stability. Hence, 
in real materials, the detailed arrangement of Skyrmions 
into lattices and their stability strongly depend on ad- 
ditional magnetic couplings such as anisotropics, applied 
and dipolar stray fields, temperature, and static disorder, 
(iii) Skyrmions, both isolated and bounded, are charac- 
terized by axisymmetrical distribution of the magnetiza- 
tion in their cores (Fig. [2] (b-f)). 

The main features of chiral Skyrmions highlighted in 
this paper provide a basis for detailed analysis of recent 
experimental findings 0, H, Q and allow to establish re- 
lations between the classical phenomenological model (1) 
and other theoretical approaches (see [2|, [5[ and bibliog- 



raphy of these papers). 
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